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ABSTRACT 

We discuss the rotational excitation of small interstellar grains and the 
resulting electric dipole radiation from spinning dust. Attention is given 
to excitation and damping of grain rotation by: collisions with neutrals; 
collisions with ions; "plasma drag"; emission of infrared radiation; emission 
of electric dipole radiation; photoelectric emission; and formation of H2 on 
the grain surface. Electrostatic "focussing" can substantially enhance the rate 
of rotational excitation of grains colliding with ions. Under some conditions, 
"plasma drag" - due to interaction of the electric dipole moment of the grain 
with the electric field produced by passing ions - dominates both rotational 
damping and rotational excitation. We introduce dimensionless functions F and 
G which allow direct comparison of the contributions of different mechanisms to 
rotational drag and excitation. 

Emissivities are estimated for dust in different phases of the interstellar 
medium, including diffuse HI clouds, warm HI, low-density photoionized gas, 
and cold molecular gas. Spinning dust grains could explain much, and perhaps 
all, of the 14 - 50 GHz background component recently observed by Kogut et 
al. (1996), de Oliveira-Costa et al. (1997) and Leitch et al. (1997). Future 
sensitive measurements of angular structure in the microwave sky brightness 
from the ground and from space should detect this emission from high-latitude 
HI clouds. It should be possible to detect rotational emission from small grains 
by ground-based pointed observations of molecular clouds. 



Subject headings: ISM: Atomic Processes, Dust, Radiation; Cosmic Microwave 
Background 
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1. Introduction 

Experiments to map the cosmic microwave background radiation have stimulated 
renewed interest in diffuse Galactic emission. Most recently, Kogut et al (f996), de 
Oliveira-Costa et al. (f997) and Leitch et al. (1997) have reported a new component 
of galactic microwave emission which is correlated with 100/xm thermal emission from 
interstellar dust. 

Kogut et al. found the emission excess to have Jj,(31.5GHz) Jj,(53GHz), and Leitch 
et al. found /„(14. 5GHz) ps /^(32GHz), consistent with the spectrum of free-free emission, 
but nondetection of Ha emission in these directions is inconsistent with free-free emission 
accounting for the microwave excess unless the plasma temperature T > 10 6 K. Leitch et al. 
therefore proposed that the observed emission was free-free emission from shock-heated gas 
in a supernova remnant. 

Draine & Lazarian (1998, hereafter DL98) showed, however, that the observed 
microwave excess could not be due to free-free emission from hot gas, as this would require 
an energy injection rate at least 2 orders of magnitude greater than that the energy input 
due to supernovae. DL98 showed that the microwave excess could in fact be electric dipole 
emission from rapidly-rotating dust grains. 

To predict the intensity of dipole emission one needs (1) the numbers of small grains, 
(2) their dipole moments, and (3) their rotational velocities. The observed intensity of 12 
and 25/zm emission from interstellar clouds allows us to estimate the numbers of very small 
grains (Leger & Puget 1984; Draine & Anderson 1985; Desert, Boulanger, & Puget 1990). 
To estimate the dipole moments, we consider the likely displacements between the charge 
and mass centroids for grains, plus the intrinsic dipole moment arising from polarized 
chemical bonds in the grain material. However, the main thrust of the present paper is a 
comprehensive study of the rotational dynamics of small grains, in order to estimate their 
rotation rates. 

Ferrara & Dettmar (1994) showed that small grains rotating with thermal rotation 
rates could produce detectable radio emission, but they did not address the details of 
rotational excitation and damping. Rotational excitation of small grains (= large molecules) 
has been discussed previously by Rouan et al (1992), who described the effects of collisions 
with gas atoms and absorption and emission of radiation. In the present work we reexamine 
this problem, and include the important effects of collisions with ions, which were neglected 
in the study by Rouan et al. We also include the effects of "plasma drag" due to interaction 
of the grain with passing ions, first considered by Anderson & Watson (1993). We derive 
rates for rotational damping and excitation due to plasma drag which are somewhat larger 
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than estimated by Anderson & Watson (1993). 

We discuss rotational excitation as a function of both grain size and environmental 
conditions. Electrostatic focussing of ions makes them very effective at delivering angular 
momentum to the grains; ion collisions can dominate the rotational excitation of very small 
grains even in predominantly neutral regions. 

For our adopted population of very small grains, we predict the microwave emissivity 
of various phases of the interstellar medium, ranging from diffuse gas to molecular clouds. 
We expect detectable levels of emission from spinning dust grains in all of these phases. 

The paper is organized as follows: After describing the environments which we consider 
(§0), we discuss the electric dipole moments expected for small dust grains (§|j). The various 
rotational damping processes are reviewed in §[|, and the rotational excitation mechanisms 
in §|5|. Using these rates for rotational excitation and damping, we calculate the resulting 
rate of grain rotation in The importance of the impulsive nature of the rotational 
excitation is discussed in §[7|, and the effects of centrifugal stresses in §|8|. 



The reader interested primarily in the predicted emission may wish to skip j 
and proceed directly to §|9], where we describe our assumptions concerning the grain size 
distribution in various environments. We present the microwave emission spectra expected 



for grains in both diffuse regions and molecular clouds in §|10|. The detection of this 



microwave emission from pointed observations of dense clouds is discussed in §11. The 



principal uncertainties in our estimates are discussed in § |12| , and our results are summarized 

in am 



2. Environments 

We will consider grains in five different idealized "phases" of the interstellar medium: 
"Cold Neutral Medium" (CNM), "Warm Neutral Medium" (WNM), "Warm Ionized 
Medium" (WIM), "Molecular Cloud" (MC), and "Dark Cloud" (DC). In Table Q we give 
the adopted values of nn (H nucleon density), T (gas kinetic temperature), % (starlight 
intensity relative to the average starlight background), and other properties of these phases 



3. Grain Properties 
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3.1. Grain Sizes and Shapes 

We will characterize the grain size by the "volume equivalent" radius a = 10~ 7 a_7cm, 
the radius of a sphere of equal volume. We will assume a density p ~ 2. gem" 3 , having 
in mind a carbonaceous material with C:H^3:1; for C:H::3:1 the mass per atom is 
m ~ 9.25 amu, and the number of atoms in the grain is then 

N = 545a 3 7 . (1) 

We will consider grains down to a minimum size N ~ 25, or a_7 0.36. 

We require the moment of inertia as a function of grain size, so assumptions must 
be made concerning the grain shape. The smallest grains may be chainlike, as favored 
by Thaddeus (1995), or sheetlike, as expected for polycyclic aromatic hydrocarbons. For 
generality, we will suppose that grains with a < a\ are cylindrical with diameter d, grains 
with a\ < a < a 2 are disklike with thickness d, and grains with a > a 2 are spherical; we 
take d = 3.35 x 1CT 8 cm, the interlayer separation in graphite. The actual geometry of small 
grains is uncertain. Here we usually take a x = and a 2 = 6 x 1CT 8 cm (grains with N < 120 
atoms are disklike), but we will also consider the possibility that oi = 4 x 10~ 8 cm (grains 
with N < 50 atoms are linear). 

The physics of rotational excitation and damping will be explicitly derived for spherical 
grains; when nonspherical grains are considered, we will attempt to use results for spherical 
grains but with the radius a for the sphere replaced by an appropriate length scale for the 
nonspherical shape. For convex shapes, the rates for collisions of neutral grains with atoms 
depends on the surface area S, and we therefore define the "surface-equivalent" radius a s , 
the radius of a sphere with equal surface area: 

Aita 2 s = S . (2) 

The excitation of rotational kinetic energy by collisions is proportional to J r 2 dS, where 
r is the distance from the surface point to the center-of-mass. We therefore define the 
"excitation equivalent" radius a x , the radius of a sphere with equal / r 2 dS: 

Anal = J r 2 dS . (3) 

Another important quantity is the ratio of the moment of inertia to the value for a sphere 
of radius a: £ = I /(O.AMa 2 ), where M is the grain mass and / is the largest eigenvalue of 
the moment of inertia tensor. 
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3.2. Grain Charge 

Grains acquire charge through photoelectric emission and collisions with electrons and 
ions. Bakes & Tielens (1994) have discussed the rate N pe (Z) of photoelectron emission 
from a carbonaceous particle with charge Ze; we adopt their rates with minor modifications 
(Weingartner & Draine 1998). Collisional charging processes have been discussed by Draine 
& Sutin (1987), whose rates we employ here. When considering nonspherical grains, we 
approximate the potential as U ~ Ze/a s . In computing the electron and ion capture rates, 
we assume a spherical geometry with a radius a s . 

Let N e (Z) and N%(Z) be the rate of electron capture and ion impact for a grain of 
radius a s and charge Ze. With these rates, we solve the usual equations to obtain the 
steady-state charge distribution function f(Z) for grains of radius a: 

N i (Z) + N pe (Z)]f(Z) = N e (Z + l)f(Z + l) . (4) 

Selected charge distribution functions f(Z) are shown in Figure [l] for CNM and WNM 
conditions. 

The charge distribution for grains of radius a has mode Z m , centroid (Z), and variance 
cr| = (Z 2 ) — (Z) 2 ; these are shown in Figure 0. The time scale for changes in the grain 
charge is of interest. We define a characteristic charging time by 

T = 1 + 4 , 5) 

TZ N e (Z m ) + Ni(Z m ) + N pe (Z m ) 

When cr| -C 1, Tz is simply the lifetime of charge state Z m ; when cr| ^> 1, Tz is an estimate 
for the time required for the grain charge to change by of order the standard deviation oz 
of the charge distribution. 

Figure |3| shows Tz as a function of grain size a for CNM and WNM conditions. 



3.3. Quantum vs. Classical Dynamics 

The grain angular momentum Iuj is quantized. If the rotational kinetic energy is 
(3/2)kT rot , then the rotational quantum number 

J = ^ = 5.85£ 1/2 iV 5 / 6 ( ^™L) . (6) 



h * Viook/ 

The smallest grain we consider will have £ 3//5 iV £ 35, so that J £ 113(T rot /100K) 1 / 2 3> 1. 
We are therefore justified in using classical mechanics to discuss the grain rotation. 
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3.4. Grain Dipole Moments 



A spinning grain radiates power 



P 



2 uj V 2 sin 2 6 

3 7 3 



(7) 



where 6 is the angle between the angular velocity u> and electric dipole moment fi. We 
assume that the grain dipole moment fi can be written 



where Hi is the intrinsic dipole moment of an uncharged grain, Ze is the grain charge, and 
the vector ea x is the displacement between the grain center of mass and the charge centroid 
[the length a x is defined by eq.([|)]. As noted by Purcell (1976), an irregular grain, even 
if perfectly conducting, will, in general, have its charge centroid displaced from the mass 
centroid. The magnitude of the displacement depends on the grain shape; we will assume it 
to be ~ 10% of a x , or |e| ~ 0.1. Thus 



Neutral molecules can of course have electric dipole moments. Except for the case 
of highly symmetric molecules, these dipole moments are expected to be appreciable. If 
the grain material has short range order, but minimal long-range order, then it consists of 
more-or-less randomly-arranged chemical substructures, and 



To estimate the likely magnitude of /3, we show in Table |2| the dipole moments associated 
with particular bonds in hydrocarbon molecules. Dipole moments of ~ 1 debye are evidently 
typical in hydrocarbon molecules. Of course, a highly symmetric molecule would have a 
net dipole moment of zero, but we do not expect such perfect symmetry to be the norm for 
large interstellar molecules. The larger molecules should have minimal long-range order, 
having grown by random arrival of gas atoms or by coagulation with other small molecules. 
It might be argued that very small molecules will be predominantly symmetric because 
only the most robust molecules will survive in the interstellar ultraviolet radiation field. 
For example, the highly symmetric polycyclic aromatic hydrocarbon C24H12 coronene (with 
Hi = 0) has been considered as a candidate for the very smallest interstellar grains. We 
note, however, that under interstellar conditions the carbon skeleton for very small grains is 
expected not to be fully hydrogenated (Omont 1986). Such radicals are not easily studied 



H = Hi + eZea x 





(9) 



I* « N 1/2 (3 



(10) 
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in the laboratory, and dipole moments are unavailable. We can estimate, however, that 
substitution of OH for one of the H in coronene would produce a molecule with a dipole 
moment of ~ 1.6debye, or j3 ~ 0.3debye. In Table [3] we give electric dipole moments for 
selected molecules, as well as the corresponding values of (3 = fii/N - 5 . 

Based on Table |3|, and the expectation that radicals will typically have larger dipole 
moments than complete molecules, we will take (3q = 0.4 debye as representative.^ 

In an ensemble of grains of radius a, we may suppose the direction of e and the 
direction of the grain's intrinsic dipole moment ^ to be uncorrelated. Thus, with e = 0.1 
and j3 = 0.4 debye the mean square dipole moment will be 

/ , 2 = [(4.8) 2 (^) 2 (Z 2 ) + (9.3) 2 a_ 7 ]a 2 _ 7 (debye) 2 . (11) 

If the orientation of n is uncorrelated with the angular velocity then (sin 2 9) = 2/3, and 
the expected radiated power per grain is just 

9 c 3 

4. Rotational Damping 

Gas-grain interactions, plasma-grain interactions, infrared emission, and radio emission 
damp the rotation of small grains. In this section we present the rates for these processes. 

We consider a spherical grain with temperature Tj in gas of temperature T, H nucleon 
density hh = n(H) + 2n(H 2 ), and with n(He) = O.Itih. 

In thermal equilibrium at temperature T, the grain would have an rms rotation 
frequency 

^1 = ~ 5.60 x 10 9 al? /2 r 1/2 r 2 1/2 Hz, (13) 

2n 

where T 2 = T/100K. If the dust temperature Td ^ T, or there are sources of rotational 
excitation and drag other than the gas (see below), then the grain rotation rate will differ 
from eq.([L3|). 

Here we summarize the contributions of gas-grain interactions to both rotational 
damping and excitation; details may be found in Appendix |B[ It is convenient to 
"normalize" the various drag processes to the drag which would be produced by "sticky" 



1 Rouan et al. (1992) assumed /i = 3 debye for N = 90, or (3 = 0.32 debye. 



(12) 
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collisions in a pure H gas of density n-^: thus for drag process j we define the dimensionless 
quantity Fj such that the contribution of process j to the drag torque is 



~dt 



( 8kT \ 2na^.m-R 



ujF, 



(14) 



\7nri11 J 

Neutral grains in a gas of pure H atoms would have F — 1, with a rotational damping time 

(15) 
(16) 



T H 



4£pa 5 /7rm H \ 1/2 
hn^muat \ 8kT ) 



3.30 x 10 10 £ 



20 cm 



-3 s 



T, 



-1/2 



n H 



0-7 — 



where we have assumed p = 2 gem 3 . We will use th as a fiducial time scale for the different 
sources of rotational damping. Thus the rotational damping time due to process j is 



The linear drag processes are additive: 



F n + F + F v + Fi 



1R 



(17) 



(1* 



where F n , Fi, F p , and F^ are the contributions from neutral impacts, ion impacts, plasma 
drag, and thermal emission of infrared photons. 



4.1. Collisional Drag 

Collisional damping is a well-studied process (see Jones & Spitzer 1967). We assume 
that when species (e.g., H, H 2 , He) arrive at the grain surface they temporarily "stick" 
and then are desorbed with a thermal velocity distribution relative to the local (moving) 
grain surface (the possible effects of H 2 formation are considered separately in §|5.6|). The 
dimensionless factor Fj for species j, which would be unity for H atoms and neutral grains, 
depends on the grain radius a, charge state Z g , and on the gas composition and temperature 
T. 

Let f(Z g ) be the probability that the grain will have a charge Z g e. Because the 
charging timescale tz is short compared to the timescale tj over which the grain rotation 
rate changes appreciably, we can average over the grain charge distribution f(Z g ). Then 
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(cf. Appendix |Eg and 



nn_ fmn\ 1/2 



z n n n Vm H / 



exp(-Z 2 e 2 n ) + \Z g \e n 7!- 1/2 eif(\Z g \e n ) 



(19) 



F; 



where 



ttr7\^ n i ( mi\ 1 / 2 ( ' ZgZie 



+ /( z^o)E^p) 1/2 (i + ^, 



1 — x if x < 
e~ x if x > 



2 6 ^"71 

6n = 2^fcT 



2e 



2 \ 1/2 



(20) 
(21) 

(22) 
(23) 
(24) 



n n , m n and a n are the number density, mass, and polarizability of neutral species n. 
Equation ( |T9| ) requires summation over neutral species n, and eq. fl2"0"|) over ion species i 
F n and F< are shown in Figs, g - for MC, CNM, WNM, and WIM conditions. 



4.2. Plasma Drag 

If the grain has an electric dipole moment fi, then the interaction of this electric dipole 
moment with passing ions couples the grain rotation to the plasma (Anderson & Watson 
1993). The contribution to the damping is derived in Appendix IB.^I 

Fp "3^F?i^(^) {MUa) + cos 2 v&ln[min(VA D /M]} , (25) 
where ^ 

° V 4nn P e 2 j ^ 



For nonspherical grains we use the length scale a s to estimate the potential ~ Ze/a s and surface electric 
field ~ Ze/a 2 s . 

3 In eq.(E3) the factor a~ 4 enters because of the normalization (Rh. 
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is the Debye shielding length, 

b„ = {I/mY' 2 , (27) 

b q = I{2kT/m l ) 1/2 /h (28) 

and \I/ is the angle between fi and u>. Our result Q25D differs from that of Anderson & 
Watson (1993), who neglect impact parameters with b > b^. In the present calculations we 
set cos 2 $ = 1/3, appropriate if ui is randomly oriented relative to //fj F p is shown in Figs. 

i-0- 



4.3. Infrared Emission 



A more subtle process is damping due to infrared emission (Purcell & Spitzer 1971; 
Martin 1972). As shown in Appendix 0, emission of a thermally-excited photon of energy hv 
from a spherical grain rotating at angular velocity u <C v removes, on average, an angular 
momentum Too j{2iiv). The grain is heated by absorption of photons from the background 
starlight, with energy density u*, which we assume to have the spectrum estimated for 
interstellar starlight by Mezger, Mathis, & Panagia (1982) and Mathis, Mezger, & Panagia 
(1983). The infrared emission cross section is assumed to scale with frequency as 

Qu = Qo(-) 2 , (29) 

as appropriate for graphite, silicate, and various other candidate grain materials (Draine 
& Lee 1984). Since there is reason to believe that the smallest interstellar grains are 
predominantly carbonaceous,^ here we will assume optical properties appropriate for small 
graphitic particles (though the discussion in Appendix |C] is fully general). 

If the grain is large enough that its temperature may be approximated as a constant 
Ta, then the rotational drag is characterized by (see Appendix |C.l|) p| 

FlR ^^f^f(^)(^£f(^' . ,30) 
a-7 Vmisrf/ V n H J \ T J \a x J 



4 We note, however, that disklike grains would tend to have cos 2 ^ w 1, since internal relaxation will 
cause the principal axis of largest moment of inertia to align with the angular momentum, whereas \x will 
tend to be perpendicular to the principal axis of largest moment of inertia. 

5 The 10/^m silicate feature does not appear in the emission spectrum of diffuse interstellar clouds (Onaka, 
ct al. 1996), implying that the particles small enough to be heated to Td XL 150K by single-photon heating 
do not have a significant silicate component. 

6 The IR drag torque depends on Qa 2 oc a 3 . The factor {a/a x ) A enters in ( p0| ) and (|3l]) from the 
normalization (|i"4|). 
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For small grains, the quantized heating by starlight photons results in thermal "spikes" 
followed by intervals of cooling (see, e.g., Draine & Anderson 1985), thereby altering the 
emission spectrum and affecting the rotational damping, as noted by Rouan et al (1992). 
As shown in Appendix (J. 2 , if individual heating events may be treated separately, for very 
small carbonaceous grains we obtain (cf. eq. C26Q 

1/2 ( u * \ /20cm" 3 \ i /2 ( a \ 4 



4.49a^ — — T 2 -^ - • (31) 
V^isrf/ \ n n J \a x J 

We have i*iR, g = -^ir, c for a = aiR, where 

/ 11 \ -2/9 

a lR = 5.68 x 1(T 7 — — cm . (32) 
Vmisrf/ 

We use eq.(^) for a > aiR, and eq. (|3T|) for a < oir, or 

Fir = min(F IRi g, F IRiC ) . (33) 
In Figs. |] - ^ we show F IR as a function of grain radius a. 



4.4. Electric Dipole Damping 



The spinning grain will radiate power as given by eq. (|12[). The associated rotational 
damping time depends on oj. It will prove convenient (see §||) to define a characteristic 
damping time 

3/ 2 c 3 16n 2 e p 2 a 10 c 3 
Ted ~ 4/i 2 fcT~ 75 fi 2 kT { ] 

= 7 ' 13 X 10l ° S ((a x /a) 2 (Z 2 ) + 3.8~(^/0.4D)2a_ 7 )T 2 ' (35) 

where we have assumed the dipole moment to be given by eq. (|TT|) . With r ed defined by 
eq. (|35"D , we see that electric dipole damping contributes 

\uodt) ed 3kTr ed ' 1 ' 

In Figures § - we show r H /r ed as a function of a. 
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4.5. Relative Importance of Damping Mechanisms 

Although all the processes discussed above contribute to rotational damping, some 
of them are more important than others. In diffuse molecular clouds (see Fig. |j) the 
plasma drag is dominant. This result is counterintuitive as the degree of ionization is small 
(x = 10~ 4 ). To understand this, recall that F p oc a _5 T~ 2 , so that for small grains and 
low temperatures plasma drag can dominate the gaseous drag, even when the fractional 
ionization is low. Dipole emission drag is of only marginal importance because of the 
relatively high gas density in the molecular region. 

In cold neutral media (see Fig. |5]) plasma drag F p dominates in the range 
6 x 10~ 8 cm ;$ a £j 1.5 X 10~ 7 cm. For smaller grains damping via dipole emission is more 
important, while for larger grains infrared emission damping dominates. The infrared 
damping Fir oc a -1 for large grains (with steady temperatures) and Fir oc a 1 / 2 for small 
ones (which undergo temperature "spikes"). 

In warm neutral media (see Fig. [(]) and in warm ionized media (see Fig. [!]) the dipole 
emission is dominant for grains <^ 1.5 x 10~ 7 cm mostly as the consequence of less efficient 
coupling between grains and the surrounding gas. For larger grains infrared damping 
dominates. The difference between the two media evidently stems from the fractional 
ionization: in ionized media interactions with ions and plasma effects are more important 
than collisions with neutrals. 



5. Rotational Excitation 

5.1. Recoil from Thermal Collisions and Evaporation 

An initially stationary grain will have its rotational kinetic energy increasing at an 
average rate 

d (1 2 \ I 8kT \ 2ixa\m,Y{kT 



dt \2 



Iu 2 =n H V— (Gn + Gi + Gp + Gm) , (37) 

/ V 7rm H J I 



where the normalized excitation rate G n is due to impacting neutrals, is due to impacting 
ions, G p is due to plasma drag, and Gir is due to recoil from infrared emission. 

Neutrals deposit angular momentum when they impact the grain, and give the grain 
an additional kick when they subsequently evaporate. These processes (see Appendix |B.2| ) 
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result in a normalized excitation rate 



Gn = E/(^)E^(?) 

Z n ^ n H VMh/ 



1/2 



exp( 



± ev 

T 



exp(-^ 9 2 e^T/T e . 



2Z^ n T/T e „ 



(3« 



We separate the ion contribution into that delivered by incoming ions, and that due to 
the evaporating neutrals resulting from the ion collisions: 

d = Gf n } + Gf v) , (39) 



G : 



(in) 



z g ^o 

+/(o)E 



11: 



f m, \ 
2n H Wh/ 



1/2 



1 + 



37rV2 j 
4 ^ 2 



(40) 



G 



(ev) 



z 9 ^o t 2n R \m H / 1 



exp( 



exp(-Z 2 e 2 ) + |Z g |e i 7rV2 er f(|Z fl |e i ) _ 



+/(o)E 



2riH Vmn/ T 



1/2 r 



1 + 



7T 



1/2 



1 - x + x 2 /2 if x < 
e~ x if x > 



(41) 



(42) 



where we once again sum over neutrals n and ions i. <p is defined by eq. (|24|) and e, is defined 
by eq. (p3|) but using the polarizability of the species resulting from neutralization of the 
incoming ion i. In Figs. p|~|TT1 we show G n and Gi as functions of grain radius a. Coulomb 
attraction makes Gi the dominant excitation process for small grain. 



5.2. Excitation by the Plasma 

Excitation of grain rotation by the fluctuating electric field from passing ions is [cf. 



G » -Ztu^ E ~~ f ^z) V2 {Hbja) + cos 2 * In [min(6 (?) X D )/b u 



3a^(kT) 2 i riu v^h 



(43) 
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5.3. Infrared Emission 

Each infrared photon carries away angular momentum h, and hence there must be a 
corresponding change in the angular momentum of the grain. The energy in the radiated 
photons is due to absorption of starlight photons, but since there are ~ 500 infrared photons 
emitted per starlight photon, we may neglect the angular momentum change due to the 
photon absorption event itself. 

If the grain temperature is approximated as constant and we assume graphitic grains, 
then (see Appendix 

GlRf= m(±Y(^f(^Li) T -^ . (44) 

0-7 \a x J Vmisrf/ \ n n J 

For very small grains, the grain temperature may be treated as zero except immediately 
following a photon absorption event, in which case for graphitic grains we find: 

^ = ^(-) 4 (^)(— )^ • («) 

We take 

G m = min(Gi RjC , G IR)? ) . (46) 



5.4. Photoelectric Emission 

The contribution to rotational excitation by photoelectrons emitted randomly from the 
grain is 

°- = ^J^alkT g W»- (<^> + ( -^) ■ (47) 

where N pe (Z g ) is the rate of emission of escaping photoelectrons from a grain with charge 
Z g e, and (E pe ) is the mean kinetic energy of the escaping photoelectrons at infinity. We 
assume the photoelectrons emerge from the grain surface with a "cosine-law" directional 
distribution. (E pe ) is as estimated by Bakes & Tielens (1994). Because m e <C m H , the 
contribution of photoelectric emission to rotational drag is negligible. 



7 The factor (a/a x ) A in eq.(|44|) and ( p5| ) enters from the normalization of G [see eq.(|37|)]. 
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5.5. Random H 2 Formation 

If a fraction 7 of arriving H atoms recombine at random on the grain surface, and 
the nascent H 2 has an average translational kinetic energy Ef, then the contribution to 
rotational excitation is 

G »*-I {1 - y) kf[ 1+ 2m H E t al ) ' (48) 

The nascent H 2 is assumed to leave the grain surface with a "cosine-law" directional 
distribution, with the internal angular momentum Jh uncorrelated with the velocity. 

The efficiency 7 of H 2 formation on very small grains is very uncertain, as is the 
effective number N r of active recombination sites. For very small grains, H 2 formation 
may be suppressed due to temperature fluctuations (which limit the residence time 
of H atoms) or chemical "poisoning" (Lazarian 1995) of recombination sites, but the 
details are very uncertain. The empirical rate coefficient for H 2 formation on grains 
R Pd 3 x 10 _17 r7,(H)riH cm 3 s _1 (Jura 1975) provides a limit: if small grains contribute a 
geometric cross section per H ~ 2 x 10~ 21 cm 2 /H (see §|9|), then 7 ;$ 0.2 if H 2 formation 
on small grains is not to exceed the empirically-determined total H 2 formation rate in the 
CNM. 

In Figs. §^TT| we show Gn 2 for 7 = 0.1, E t = 0.2 eV , and (J(J + 1)) = 10 2 . We see that 
Gh 2 does not make a major contribution to rotational excitation. 



5.6. Systematic Torques 

Superthermal rotation may result if systematic torques act on a grain. Purcell 
first showed that various processes, including variations of the accomodation coefficient, 
photo-electric yield, and formation of hydrogen on preferential sites, can drive suprathermal 
rotation. Recently, radiative torques were identified as an important means of suprathermal 
spin-up (Draine & Weingartner 1996). 

Thus far, suprathermal rotation has been discussed mostly in relation to grain 
alignment and for relatively large grains, i.e. with a > 10~ 6 cm (see Lazarian & Draine 
1997). The degree of suprathermality, i.e. the ratio of grain kinetic energy to its energy of 
thermal rotation, depends on grain size. 

Systematic torques differ from the random torques in having a nonzero time average 
(in grain coordinates). Let T s be a time-averaged systematic torque in grain coordinates. 
If the grain were not subject to fluctuating torques, it would attain a steady rotational 
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velocity u s satisfying 



2, ,3 

s 



th 9c 3 

For our purposes it is sufficient to use the simple approximation 



IF 



1 + 



v 

9r,c 3 



r s r H \ 3 

IF J 



-1/3 



(49) 



(50) 



which is exact in the limits r —> and T s 
does work on the grain is then just T s u s . 



oo. The rate at which the systematic torque 



The actual grain rotational velocity u is also the result of fluctuating torques, but T s u s 
remains a good estimate for the mean rate at which the systematic torque contributes to the 
rotational kinetic energy. Thus it is natural to characterize the importance of superthermal 
torques through a dimensionless quantity G s defined by analogy with the dimensionless 
excitation rates G n , Gi, and G p : 



G s 



n u (2inn H kTy/Ha*kT 



(51) 



where u s is given by eq. (|50|). 



Formation of H2 molecules on catalytic sites can drive superthermal rotation. Averaged 
over grains of radius a, the rms torque is0 



4-fn(H.)a 2 x a s 



N 



1/2 



(Trm-iiEfkT) 



1/2 



(52) 



where 7 is the fraction of impinging H atoms which leave the grain surface as H2, 
Ef = 0.2i?feV is the translational kinetic energy per nascent H2, N r is the number of 
recombination sites on the grain surface, and we have assumed that at each recombination 
site the newly-formed H 2 leaves the surface in a fixed direction, with the direction for each 
site randomly-drawn from a "cosine-law" distribution. Then, using eq. (|50| , ^Tl) 



G, w 1.75 x 1(T 



N r 



FTo 



vf 
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(53) 



For large grains, we expect N r oc a 2 , but for all grains where H 2 forms (7 > 0) we must 
have N~ > 1. We will consider two extreme cases: 



The rms moment-arm oc (a^/a^) 1 / 2 , and the H2 formation rate per site oc a 2 s /N r , hence T s oc a^a s . 
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1. No H 2 formation: 7 = 0, T s = 0, and G s = 0. 

2. H2 formation with 7 = 0.1, E{ = 1, and the number of recombination sites on a grain 
taken to be 

iV r = 1 + ( % ) . (54) 

V2xl0- 7 cm/ v ; 



5.7. Relative Importance of Excitation Mechanisms 

Unlike the case of the damping mechanisms, ion bombardment dominates rotational 
excitation for small grains in all ISM phases discussed above. This apparent disparity 
between damping and excitation processes, which seemingly contradicts the Fluctuation- 
Dissipation theorem, is the consequence of the non-equilibrium state of the interstellar 
medium. The ionization in T £ 10 4 K regions is primarily due to photoionization and cosmic 
rays, and far exceeds the ionization which would correspond to thermodynamic equilibrium 
at the kinetic temperature of the gas. The dominance of ion excitations stems from the 
increase of both capture cross section and angular momentum per capture. 

In molecular clouds (see Fig. ||) excitation through ion collisions is dominant for grains 
^ 1.7 x 1CT 7 cm. As a result, these grains rotate with velocities a bit larger than would 
result from Brownian motion. 

In cold neutral media the picture is more involved (see Fig. |S|). Regular torques 
associated with H2 formation can affect grain dynamics for a ^ 10~ 7 cm. The rate of 
excitation through infrared emission is very close to that of gaseous bombardment in the 
range 4 x 10 -8 cm< a < 4 x 10 -7 cm and only for larger grains do grain-neutral interactions 
dominate grain excitation. 

The relative importance of ion excitations drops with grain temperature. Therefore 
for warm neutral media (see Fig. ITUI) collisions with ions constitute the leading excitation 



mechanism only for grains less than ~ 10 6 cm. For larger grains collisions with neutrals 
dominate. However, for warm ionized media (see Fig. |TT]) , where 99% of species are ionized, 



collisions with ions are bound to be the major source of rotational excitation. 
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6. Rate of Rotation 

With F defined by equation flI3|), the time-averaged rate of change of rotational kinetic 
energy may be written 

\dt 2 I r H r H 9c 3 x 1 ' v ; 

where 

G = G n + G 4 + G p + G IR + G s , (56) 

where G s > if there are systematic torques due to H2 formation [eq. (|53|) 1. Although the 
distribution of angular velocities will not be Maxwellian (except in the limit where electric 
dipole emission and superthermal torques are both negligible), we assume the relation 
between moments of a Maxwellian distribution, 

{^) = \{^) 2 , (57) 



to obtain a quadratic equation in (u> 2 ) from eq. flSBp. With the electric dipole damping time 
r e d defined by eq. fl35|), the solution is 



{Uj2) ~ 1 + [1 + (G/^)(20r H /3r e(i )]V 2 (f) (T) " (58) 

Eq. fl58|) provides a reasonable estimate for the rms grain rotation rate. 

The effective rotational damping time tj is given by 

1 F /iw 2 \ 1 , N 

- = -+-—— • (59) 
tj r H \3kT J r ed 



7. Effects of Impulsive Torques 

The above discussion of grain rotation has been directed at estimating (u; 4 ) for grains 
of radius a. Because we do not have thermal equilibrium, however, the angular velocity 
distribution may depart substantially from a Maxwellian. In particular, if I5u is the 
angular momentum of an impacting particle, then if (Su 2 ) £ (uj 2 ), then the angular velocity 
distribution will depart from a Maxwellian. This occurs when infrared emission and/or 
electric dipole emission is able to reduce the grain angular velocity substantially between 
impact events. 
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One can show that colliding species j has 

({Suj^j = (G J /F j )2a 2 m J kT/I 2 . (60) 

In Figure O we show {(Su) 2 ) for impacting ions, normalized by (u 2 ) for the grain. We 
see that for large grains, the impulsive nature of the rotational excitation is unimportant, 
and the time evolution of the grain rotation can be treated as a continuous process. For 
grains with a £ 7 A, however, a single ion impact can change the grain angular momentum 
substantially. From Fig. || we see that ion impacts are the dominant source of rotational 
excitation for small grains in the MC, CNM, and WIM phases. This is primarily the result 
of the strong electric dipole damping for small grains (see Figs. || - 0) in diffuse gas, which 
causes the smaller grains to have sub-thermal rates of rotation. 



8. Centrifugal Stresses 

An upper limit for the grain rotational frequency u can be obtained from the ability of 
grains to withstand centrifugal stress (Draine & Salpeter 1979). With a characteristic stress 
(l/4)puj 2 a 2 x) and a maximum stress Sjyiqxi the maximum rotational frequency is 




(61) 

Bulk polycrystalline substances have S max ~ 10 9 ergs cm -3 , while for single- crystal materials 
Smax ~ 10 11 ergs cm" 3 . As can be seen from Figure O, the expected rotation rates are such 



that even grains with S max as small as 10 ergs cm can survive 



9. Size Distribution 

9.1. Diffuse Clouds 

As shown below, microwave emission from grains is only significant for uo/2-k £ 1 GHz, 
so from Fig. |l| we see that only a^3x 10~ 7 cm grains (N £ 10 4 ) are of capable of 
contributing to such emission. 

The size distribution of these ultrasmall dust grains is very poorly known. Studies of 
interstellar extinction are relatively insensitive to the detailed distribution of very small 
dust grains. Observations of extinction at A £ 1000 A constrain the total mass in dust at 
fl^lx 10~ 6 cm, but do not reveal how this mass is distributed over grains of different sizes. 
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The MRN distribution (Mathis, Rumpl, & Nordsieck 1977), extended down to very 
small sizes, provides one estimate for the population of very small grains. The graphitic 
component has 

dn 

— = U}iA MRN a~ 3 ' 5 for a min < a < 2.5 x 10~ 5 



da 



cm, 



(62) 



with A 



MRN 



10 



-25.16 



cm 2 ' 5 (Draine & Lee 1984). This size distribution, if extended down 
to a min = 3.6 A, has 17% of the grain mass in a < 10~ 6 cm grains, but only 2.6% of the 
mass is in a < 10~ 7 cm grains. This distribution of graphitic grains would contribute a 
geometric cross section per H atom Emrn = 2.2 x 10~ 21 cm 2 . 

Observations of the "UIR" emission features, as well as the strong 12/im and 25/rni 
emission observed by IRAS (Boulanger & Perault 1988), have been interpreted as showing 
that there must be a substantial population of very small grains (Leger & Puget 1984; 
Draine & Anderson 1985). Here we consider a log- normal size distribution for this 
population, which we add to the MRN distribution: 



1 dn _ 35 

— -j- = Am rno- 
rin da 



Ba 1 exp 




(63) 



We will consider a = 6 A and a = 0.4; The coefficient B is chosen so that the grains 
contain a fraction fc = 0.05 of the total carbon abundance (C/H = 4.0 x 10~ 4 : Grevesse et 
al. 1991). This grain component then contributes an additional geometric cross section per 
H atom E = 1.7 x 10~ 21 cm 2 . 

For comparison, we note that the grain model of Desert, Boulanger, and Puget (1990) 
has fc = 0.09 in the "PAH" component, molecules with 60 < iV < 540 C atoms. 



9.2. Dense Clouds 

The abundance of ultrasmall grains in dense interstellar clouds is not well known. In 
denser regions, the wavelength X m ax of maximum polarization is increased, and the ratio 
Ry = Ay/E(B — V) ~ 5.5, significantly larger than the value Ry ~ 3.1 characteristic of 
diffuse regions; both are indications that the fraction of the grain mass in smaller grains is 
reduced in denser regions. Here we will assume that the fraction of the mass in ultrasmall 
grains is reduced by a factor 5 relative to diffuse clouds: for a £ 3 x 10~ 7 cm we adopt the 
size distribution flB3]) but with A MRN and B reduced by a factor 5 from the values used for 
diffuse clouds. 
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10. Emissivity 

Recognizing that there will be a range of dipole moments among grains of a given size, 
we will assume that 25% of the grains have (3 = 0.5/3 , 50% have (3 = /3 , and 25% have 
(3 = 2/3 , with Pq = 0.4debye as our standard value. 

For each grain size a and value of (3, there is a mean square angular velocity (uj 2 ). If 
we assume u to follow a Boltzmann distribution, then the emissivity per H is 



V 2 1 f dn /A; 6 /-3^ 2 \ 

rfa -J-T-^7?exp — — . (64) 



nn \37r/ riftC 3 J da (u> 2 ) 3 / 2 \2(ui 2 ) 



We evaluate this emissivity for the size distribution of Fig. 13, for conditions characteristic 
of different phases of the interstellar medium. We consider the case where no H2 formation 
takes place, as well as the case where it occurs with 10% efficiency (7 = 0.1). 

Ferrara & Dettmar (1994) predicted the microwave emission from grains under WIM 
conditions, obtaining an emissivity j u /nn = 1.1 x 10~ 15 (z//100 GHz) 2 8 Jy sr _1 . The slope 
and magnitude of this emission is similar to our results near 10 GHz, but this is coincidental: 
Ferrara & Dettmar assumed grains to be rotating with thermal rotation rates corresponding 
to ~ 10 4 K, whereas we find very strong damping by both microwave emission and infrared 
emission (see Fig. [7]) so that the grains actually rotate very subthermally (relative to the 
8000K gas temperature in the WIM). For example, for a = 10 _7 cm grains under WIM 
conditions, we estimate uj ra d/27r « 8.5 GHz (see Fig. 0), whereas Ferrara & Dettmar 
estimate a rotation rate ~ 50 GHz for this size; since the emitted power oc a; 4 , this is a 
significant difference. In Ferrara & Dettmar's model, the 10 GHz emission comes from 
grains with a»2x 10~ 7 cm, whereas we find here that the 10 GHz emission comes primarily 
from grains with a ~ 9 x 10~ 8 cm; these two sizes differ by a factor of ~ 10 in mass, and 
~ 50 in moment of inertia. 

In Figure [16] we show the effects of decreasing the intrinsic dipole moments /ij by a 
factor of two, by taking /3 = 0.2debye and e = 0.05. In Figure |17] we show the emissivity 
calculated for H 2 formation efficiencies 7 = and 0.1; the spectra are essentially unchanged. 
It is apparent that the emission spectra are not especially sensitive to the precise values of 
A) or 7. 

It is evident that the population of small grains which has been inferred previously 
from observations of both the UIR emission features and the 12 and 25/im IRAS emission 
should produce substantial electric dipole emission in the 10-100 GHz region. 

It appears that this emission has already been detected by experiments designed to 
measure angular structure in the cosmic background radiation. Kogut et al. (1996), de 
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Oliveira-Costa et al (1997), and Leitch et al. (1997) have reported detection of 14 — 90 GHz 
microwave emission which is correlated with 100/im emission, and therefore apparently 
originates in interstellar gas or dust. Draine & Lazarian (1998) argue that this emission is 
in fact due to rotational emission from spinning dust grains in diffuse interstellar gas at high 
galactic latitudes. The observed emissivity per H nucleon is shown in Figs. [15| - [16], and it 
is seen that the observed emissivities are in fact approximately equal to the electric dipole 
emission predicted here. Future measurements of the spectrum of this diffuse emission, 
and its correlation with interstellar gas and dust, will test this interpretation, although the 
emission is weak because of the relatively small amounts of dust at high galactic latitudes. 
The MAP mission, to be launched in 2000, is expected to obtain accurate maps of the 
emission from interstellar dust in five bands, from 22 - 90 GHz. 



11. Detecting Dark Clouds 

It may also be possible to detect emission from spinning dust grains in dense clouds, 
where the larger column densities result in increased intensities. As an example, we consider 
the central surface brightness for the model which Lazarian, Goodman & Myers (1997) 
adopt for the L1755 dark cloud, with a central column density Nn = 2 x 10 22 cm -2 in 
material with densities % ~ 10 4 cm -3 , plus Nn = 1.3 x 10 21 cm -2 in diffuse molecular gas 
with nn w 300 cm~ 3 . To this we add Nn = 1 x 10 21 cm -2 in HI surrounding the cloud, with 
conditions characteristic of the CNM phase. 

We add vibrational ("thermal") emission to the rotational emission, using dust 
temperatures of 12K, 15K, and 18K for dark regions, diffuse molecular gas, and CNM, 
respectively; the assumed dust opacity is oc u 1 - 7 and reproduces the observed A £ 100/im 
emission from the interstellar medium for a grain temperature Td Pd 18K. Figure [IB] shows 
the resulting spectrum. We see that at 10 — 30 GHz antenna temperatures of ~ lmK are 
expected toward the central regions; for the adopted cloud model about half of this emission 
originates in the extended "CNM" envelope of the cloud. 



12. Discussion 

We have presented above a detailed study of rotational excitation and damping for 
grains with sizes less than ~ 10~ 6 cm. In addition to processes included in earlier work 
by Rouan et al. (1992), we have included both direct collisions with ions and "plasma 
drag". Ion collisions and plasma drag, together with damping by infrared and microwave 
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emission, dominate rotational excitation and damping for most interstellar environments. 
It is therefore important to include these processes in other studies of rotational excitation, 
such as those proposing to explain features of the diffuse interstellar bands (Rouan, Leger, 
& Coupanec 1997). 

For the assumed distribution of grain sizes and electric dipole moments, we predict 
microwave emission from spinning dust grains which can account for the "anomalous 
emission" observed recently in the 15-90 GHz range. The predicted intensities are uncertain, 
however, as they depend on three poorly-known factors: 

1. The abundances and size distribution of the small grains. This uncertainty is greatest 
in the case of dense/dark regions, where the absence of starlight denies us the evidence 
(12 - 25/im emission) which requires an abundance of ultrasmall grains in diffuse 
regions. 

2. The charge distribution of the small grains, which affects the rates of rotational 
excitation and damping. We have used standard estimates for photoelectric cross 
sections (Bakes & Tielens 1994) and electron capture cross sections (Draine & 
Sutin 1987) for very small grains. The uncertainties in these quantities affect both 
the electric dipole moment and, more importantly, the rate of angular momentum 
exchange with ions. 

3. The electric dipole moments of both neutral and charged small grains. 

The shape of the small grains (chainlike vs. sheetlike vs. quasispherical) is also uncertain, 
but is less critical than the above factors. A fifth factor - the efficiency of H2 formation 
on small grains - is not critical for these estimates: in Fig. |T2| only the grains in the CNM 
show appreciable sensitivity to whether or not H 2 formation takes place, and even for the 
CNM component the emissivity at v ^ 2 GHz is only slightly changed (see Fig. IPTf) . 

Because of these uncertainties, definitive predictions for the rotational emission 
spectrum are not yet possible. Nevertheless, within the existing theoretical and 
observational uncertainties it appears that much or all of the observed 15-90 GHz 
"anomalous" emission is due to spinning dust grains. The largest discrepancy between 
observation and theory is at 14.5 GHz, where Leitch et al. report emission about 3.5 times 
stronger than the emission predicted in Fig. |T5[ Additional measurements at v £ 20 GHz 
will be of great value to clarify whether this emission has another origin, or whether some 
of our assumptions concerning the dust must be modified. 

Emission from rotating grains must be allowed for in studies of the cosmic microwave 
background radiation (see DL98). It appears that the small rotating grains may be partially 
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aligned with the local magnetic field (Lazarian & Draine 1998), so that the electric dipole 
radiation will be linearly polarized; this may present a problem for interpretation of CMB 
polarization measurements by the MAP mission. 

13. Summary 

The principal results of this paper are as follows: 

1. Even neutral dust grains are expected to usually have electric dipole moments arising 
from polarized chemical bonds within the grain. Charged grains have an additional 
contribution to the dipole moment due to displacement of the charge centroid from 
the mass centroid. Our estimate for the dipole moment is given by eq. (|TT|). 

2. The excitation and damping of rotation in small grains is determined by collisions 
with ions and neutrals, "plasma drag" , emission of infrared and microwave radiation, 
and formation of H2 on the grain surface. Ion collisions and plasma drag, omitted in 
previous estimates of rotation rates, are included in the present analysis and found 
to often dominate rotational excitation and damping. Induced-dipole attraction of 
neutrals by charged grains, and of ions by neutral grains, can also be significant. 
Because the charge state of the grain, and the fractional ionization of the gas, do 
not reflect thermodynamic equilibrium, the fluctuation-dissipation theorem does not 
directly apply. 

3. For very small grains (a ^ 7 x 10 _8 cm), the angular momentum of colliding ions 
is large compared to the r.m.s. angular momentum of the grain, and therefore the 
grain rotation history consists of "rotational spikes" separated by intervals of gradual 
rotational damping. 

4. The estimated grain rotation rates are such that the small grains which have been 
postulated to explain the near-infrared emission feature can account for the emission 
observed at 30 - 50 GHz by Kogut et al., de Oliveira-Costa et al., and Leitch et al. 
(1997). 

5. The emission observed at 14.5 GHz by Leitch et al. is stronger than we estimate for 
spinning grains by a factor ~ 4. Additional determinations of emission from dust at 
frequencies ^ 20 GHz will be of great value. 

6. We predict that dark clouds should produce detectable microwave emission, with 
antenna temperatures of ^ lmK at ~ 10 — 30 GHz. 
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A. Geometric Factors 



For cylindrical grains (a < a\) of diameter d and length 2b we have 

, 8a 3 
h= M 2 ' 



bd d 2 ] 1/2 
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27 Vd, 

where the "surface-equivalent" radius a s and "excitation-equivalent" radis a x are defined in 
eq.(||,||), and £ = J/(0.4Ma 2 ). For disklike grains (eti < a < a 2 ) of thickness d and radius b 
we have 

- ^a 3 \ 1/2 
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For spheres (a > a 2 ), a s = a x = a and £ = 1. Figure [L9| shows how a s /a, a x /a and £ vary 
with a. 



B. Rotational Damping and Excitation of Grains in a Partially Ionized Gas 

The angular momentum of a spherical grain in a partially-ionized gas changes due to 
several distinct scattering processes: 



1. impact of neutral atoms and molecules on the grain surface; 

2. impact of ions on the grain surface; 
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3. electromagnetic interaction of the electric dipole moment of the grain with passing 
ions. 

Each of these processes contributes to both rotational excitation and rotational damping. 

The microphysics of inelastic scattering of particles which impact directly on the 
grain surface is obviously complex. We idealize the problem by imagining that species 
which arrive at a point on the grain surface are re-emitted from the same point on the 
grain surface with a thermal distribution of velocities in the frame of reference of the local 
(moving) grain surface at the instant of emission. Thus we neglect possible effects of large 
centripetal accelerations of the grain surface when the grain is rotating rapidly. We further 
assume that all arriving neutrals and ions depart the grain surface as neutrals. 

With this idealization, we can now calculate the rate of rotational damping and 
rotational excitation. The drag and rotational excitation from different components is 
additive. 



B.l. General Considerations 

Consider a spherical grain of radius a and moment of inertia / = (87r/15)pa 5 , with 
charge Z g e, rotating with angular velocity uo. Atoms and molecules arrive at the surface 
with (on average) zero net angular momentum. Angular momentum is carried away by the 
gas atoms and molecules which are re-emitted from the surface, so the rate of damping is 
directly proportional to the rate at which mass arrives at the grain surface. In order to 
more easily compare different contributions, we will "normalize" to the rate of damping 
which would be contributed by only the H atoms in neutral atomic gas: 



-1 d T 

-Iuj 

u at 



8kT' 



1/2 



2a 2 



na m H - 



(Bl) 



with F — 1 for neutral grains and H atoms in atomic gas. The drag processes are additive, 
so that 

F = E F i > ( B2 ) 



where Fj is the contribution to F from process j. 

Suppose that the interaction potential is such that particles of species x with velocity 
v at infinity collide with the grain surface if and only if the impact parameter b < b max {v). 



-28- 



Then collisions with species x contribute 
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We also consider the rate at which collisions act to increase J 2 for a stationary grain, 
and again normalize to the rate which would be appropriate for H atoms in neutral atomic 
gas if the atoms left the grain as H atoms at temperature T: 
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For a given grain, we will consider the contributions Gj which various processes j make to 
the overall excitation rate G. The incoming particles contribute 
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Each arriving particle, after delivering its incoming angular momentum to the grain, 
later "evaporates" and escapes from the grain surface. We assume that evaporating 
particles have a velocity distribution appropriate for a grain temperature T ev . Because 
the evaporating particles interact with the grain (some evaporating particles will not have 
"escape velocity" and will return to the grain surface), we use the fact that the angular 
momentum distribution of escaping particles must equal that for the same species being 
captured from a gas at temperature T ev . Thus one can show that the evaporating particles 
contribute 
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where f x ev \v) is given by eq. (|B4|) but with T replaced by T ev , and b^ x {v) is the value of 
bm ax appropriate for particles with the physical properties of the evaporating species (i.e., 
neutrals even when the impinging species x is an ion). 



B.2. Neutral-Grain Collisions 

Prior to impact, the grain-neutral interaction is approximated by an induced dipole 
potential, 

1 Z 2 e 2 

U(r) = -^n^r- (B8) 
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where a n is the polarizability of the neutral atom or molecule, and Z g e is the charge on the 
grain. The trajectories of particles in an r -4 potential are well-known (see, e.g., Wannier 
1953, Osterbrock 1961). For an atom with initial velocity v, impact parameters 



b < b (v) 



' 4Zya n \ 1/4 
mv 2 I 



(B9) 



have "spiral" trajectories which pass through the origin, while for b > bo trajectories are 
"hyperbolic", with a finite distance of closest approach. A hyperbolic trajectory with b = &o 
has a distance of closest approach r min = bo/y/2. Thus we may use the condition a = b /y/2 
to define a velocity 
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Thus the contribution to F is 
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The contribution of neutral species n to G is 
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B.3. Ion-Grain Collisions 

B.3.1. Charged Gram: Z g ^0 
For a Coulomb law potential Z g Zie 2 /r, we have (Spitzer 1941) 

( for < 

W«) = | , _ ^y/ 2 ^ ^ (B15) 

^. \ mair / 2 a 

For ions and charged grains (Z g Zi ^ 0) we neglect the modifications to b max which result 
when the "image charge" contribution to the interaction is taken into consideration (see 
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Draine & Sutin 1987, and § B.3.2| below). The contribution of colliding ions to F is then 
(Spitzer 1941) 
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The contribution to G due to the arriving ions is (Anderson & Watson 1993) 
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(B17) 
(B18) 

(B19) 
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We assume the ions to depart from the grain as neutrals. Per departing neutral, we should 
have the same contribution to G as for arriving neutrals if the gas temperature were T ev 
(see eq. |B12j and |B14j ) . Thus departing neutrals contribute 



G 



(ev) 



T 

l 2T 



exp(-Ztf) + 2Z]t % 



where 



exp(-ZV 2 ) + \Z g \^l^\Z g \ti 



e 2 ai 



2a 4 kT P . 



(B21) 



(B22) 



where a« is the polarizability of the neutral obtained from ion i. 



B.3.2. Neutral Grain: Z g = 

In the case of a neutral grain, the collision rate is modified by the polarization of the 
grain by the electric field of the ion. If the grain is approximated as a perfect conductor, 
the interaction potential is 

Z 2 e 2 a 3 
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The collision rate for this interaction potential has been discussed by Draine & Sutin (1987). 
Using their eq. (B6)|] with v = 0, gives a critical impact parameter 



bmax{v) = a 



'4Z 2 e 2 \ 1/2 ' 



mv^a 



1/2 



(B24) 



so that 



n H \m R J 



: \ V2 



7T 



1/2 



The arriving ions contribute 



(B25) 
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37T 1 / 2 
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where 



2^2 



2Zfe 
akT 



and the departing neutrals contribute an additional 



(B27) 



G 4 M = — AF . 
2T 



(B28) 



B.4. Plasma Drag 



Consider a stationary grain with electric dipole moment /x. The electric field E at the 
grain due to nearby ions exerts a torque /jxE. To estimate the effects of such torques, we 
assume that an ion with impact parameter b passes on a straight-line trajectory at constant 
velocity v. To define the orientation of fi, let the direction of v define a polar axis, with 
9 the angle between v and fi, and let the azimuthal angle = when fi is in the plane 
containing both the ion trajectory and the grain. It is then easy to show that the angular 
momentum AJ exchanged between the ion and the grain has 




(B29) 



If we now average over random orientation of /x ((sin 2 #) = 2/3, (cos 2 6) = 1/3, 
(sin 2 0) = 1/2) and then over a thermal distribution of velocities v, we find (Anderson & 



9 Note the typographical error in eq. (B6) of Draine & Sutin (1987), which should read (2ex — v)(x 2 — 
l) 2 -x = 0. 
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Watson 1993) 

p 2 - ^r^^/^r^K^)' <b3o) 
= ^,(^) i/2 & ln(V6l) . 

For the lower cutoff we simply take the grain radius b\ = a. The upper cutoff is more 
problematic. It is clear that it cannot exceed the Debye length, 

/ kT \ 1/2 / T/100K \ 1/2 , n , 

since at larger distances the ion electric field is screened. However, for a grain rotating 
at angular velocity to > 0, an ion with velocity Vth = (2kT '/mi) 1 / 2 and impact parameter 
b £j = v t h will be almost unaffected by the rotating component of the grain dipole 
moment, as the torque averaged over the grain rotation will tend to zero. For a grain with 
a thermal rotation rate to ~ {2kT / 1) 1 / 2 we have 

K - ( 8 4^] = 4-5 x KT V-^ 1/2 V2 cm . (B32) 

Let \l/ be the angle between the dipole moment n and the rotation velocity us. In eq.( p31|) 
we would then replace 

/i 2 ln(6 2 /6i) -> fi 2 [ln(6 w /6i) + cos 2 * \n(b 2 /b w )] (B33) 

We differ here from the treatment of Anderson & Watson (1993) who neglected the 
contribution of impact parameters b > b u . The rotational excitation is quantized, and we 
expect the above classical estimate to fail when the characteristic frequency ~ v t h/b of the 
time- varying electric field from the passing ion varies falls below the frequency h/I of the 
J = — > 1 rotational transition, resulting in a "quantum" cutoff b q ~ Iv t h/fr, or 

b q = 4.1 x l(T 3 £a 5 _ 7 f r / 100K y /2 cm . ( B34) 
V mi/m H J 



Passing ions then contribute^] 

dJ 2 16 /2vrm,\ 1 / 2 



Hi 



dt 3 V kT 



Z 2 e 2 fi 2 {ln(b w /a s ) + cos 2 * In [min(6„ A^/M} , (B35) 



10 It is instructive to compare B35| with quantum-mechanical calculations. For CN (p = 1.45 debye) 
and T = 10 4 K the rate coefficient for J/h = — > 1 excitation by protons is (cru)o^i ~ 4 x 10 _6 cm 3 s _1 
(Thaddeus 1972), whereas taking cos 2 VP = 1 we estimate (2mh 2 )~ 1 dJ 2 /dt ps 1.5 x 10~ 5 cm 3 s _1 , too large 
by a factor ~ 4. Such agreement even for a case where quantum effects are dominant is reassuring. 



-33- 



Gp = ^W ^4^^M 1 ^/^) +cos2lI/ln [ min ^' A ^/^} • (B36) 

Having estimated the excitation rate G p , we may now invoke the fluctuation-dissipation 
theorem to see that the normalized contribution to the damping rate must be simply 
F p = G p , so that if torques from passing ions were the only torques present, the grain would 
attain a thermal rotational distribution. 



C. Damping and Excitation of Rotation by IR Emission 

C.l. Steady Emission 

Suppose that the grain at rest contains six rotating electric dipoles, each with electric 
dipole moment p and with angular velocity ui in grain body coordinates. There is one dipole 
rotating clockwise and one counterclockwise around each of the x—, y—, and z— axes. The 
dipoles emit incoherently, so that at rest the power radiated is just 



dE „ 2u 4 p 2 

— = 6 x 

dt 3c 3 



(CI) 



the emission is isotropic and unpolarized. Each dipole emits radiation which, in directions 
parallel to the dipole's rotation axis, is 100% circularly polarized. For an individual dipole 
the angular momentum loss dL^ /dt is just 



djW 2oj 3 p 2 



dt 



3c 3 



(C2) 



A system of two oppositely rotating dipoles radiates a net angular momentum if the 
whole system rotates with frequency u T along the dipoles' rotation axis: in the rest frame 
the frequencies of emission for the two dipoles are, respectively, u + u T and u — u r . 

Let the system of six dipoles rotate around the a>axis with frequency uj t <^ uj. The rate 
of loss of angular momentum associated with the single emission frequency u ^> u r is then 



dJ 
~dt 



2(u + u r ) 3 
3 



2 (u — u r 

3 ~ 



P uj r dE 

P 2 17 

uj 1 dt 



(C3) 



Using the Planck expression for thermal emissivity we find 



dJ 2 f°° Qv B v 

— = -u r a / 

at Jo v l 



dv 



(C4) 
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Qu = Qo[-Y , (C6) 



while 

A TP. roc 

= -4ttV / Q v B v du . (C5) 
at Jo 

Now suppose 

Q..=Or,( 

where we will assume (3 = 2 for interstellar dust (Draine & Lee 1984), and set 

dE 

— = -na 2 (Q)*u*c , (C7) 
at 

where w* is the energy density of starlight, and (Q)# is the grain absorption efficiency 
averaged over the starlight spectrum. Then 



(C8) 



(C9) 



dJ h 2 T(P + 2)(([3 + 2) (Q^oVc 
dt~ 47rr(/? + 4)C(/3 + 4) (kT d ) 2 

where T and ( are the usual gamma function and Riemann ^-function. Thus 

3 r(/? + 2)C(/? + 2) (Q)*u*ch 2 
m ' c (87r)y 2 T(p + 4)C(p + 4)n H (m H kTy/ 2 (akT d ) 2 

= mfj^)(^)n^f(^)' . (cl o) 

The grain temperature T rf and starlight intensity u* are related through 

he 



87r/icQ A(?r(/3 + 4)C(/3 + 4)_ 



(Cll) 



where Ao = c/vq. We consider three grain materials: graphite, "astronomical silicate", and 
a— SiC (Draine & Lee 1984; Laor & Draine 1993), for each of which the infrared emissivity 
has (3 = 2, so that T d oc uV 6 . 

The recoil from photon emission is a source of rotational excitation. For a nonrotating 
grain with radius a <C hc/kT d , 

d j2= d_N^ n2 

dt dt y ' 



where the photon emission rate is 

dN ph _ T(f3 + 3) (((3 + 3) 1 



dt T((3 + 4) C(P + 4) kT d 

Thus 



7ra 2 (Q)*w*c . (C13) 



(gWgc r(/? + 3)C(/? + 3) 

IR ' C 16n H m|/ 2 (27rA;T)WA;r d r(/3 + 4)C(/3 + 4) " 1 J 
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C.2. Thermal Spikes 

Grains undergo a sudden temperature rise following each photon absorption. For 
very small grains, the temperature history can be regarded as a sequence of independent 
"thermal spikes", separated by intervals where the grain is very cold. Rouan et al. (1992) 
discussed the rotational damping for a vibrational density of states chosen to approximate 
a PAH molecule. Here we adopt the Debye model for the heat capacity of the 3N — 6 
vibrational modes (Kittel 1972). If G is the Debye temperature, the thermal energy content 
for T < Q is 

Q—-4 / 'T 1 \ 3 

E = —(N-2)kT(-) . (C15) 



5 v ' ve, 

The absorbed energy of a UV photon first heats the grain and then is emitted in the 
infrared. From ( |C5|) and (|C8|) we obtain 

dJ _ u r h 2 r(/? + 2)CQ3 + 2) 

dE 47r 2 (fcT) 2 r(/3 + 4)C(/3 + 4) ' 1 ' 



Using (|C15 ) to eliminate T in favor of E, we integrate to obtain the angular momentum 



loss in cooling from initial energy E to T = 0: 

SJ = ~Au r E 1 ' 2 , (C17) 
h 2 /3y/ 2 r(/3 + 2)C(/3 + 2) (jV-2)V 2 

2 \5J r(/3 + 4)C(/3 + 4) (A;6) 3 / 2 ' 1 ' 

The energy E is provided by the interstellar radiation field. Therefore 

d 4- = -Alu t (°° dv ( ^) Q v {hvf 2 na 2 . (C19) 
at Jo \hv ) 



Define 



T = dv C ^ , (C20) 



o hv 



(E) T = — I dv cu v , (C21) 



M^hrH&O" ' (C22) 

1 r°° 

(Q)fe = ^ / dv cu u Q v , (C23) 
T Jo 

^^TWJ^^^)^ ■ <C24) 
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Then 



Thus 



IR, 9 



3 /3tt\ ^ 2 + 2)C(/3 + 2) h 2 {N - 2) 1 / 2 T(Q) r (E l / 2 )r Q 



8 VlOy r(/3 + 4)C(/3 + 4) a 2 (fce) 3 / 2 nnlmHfcT) 1 ^ 



(C25) 



(C26) 



Because oc a, F IR)9 oc (N — 2) 1/2 /a. 

The emission following each heating event also contributes to the rotational excitation 
of the grain, with 



AJ 2 



h 2 T(J3 + 3)t(j3 + 3) 



Thus 



G 



IR,q 



3nT([3 + 4)((f3 + 4) 

h 2 T{(3 + 3)({(3 + 3) (iV-2) 1 / 4 
37r(/? + 4)C(/? + 4) (A;0) 3 / 4 



< w - 2)1/4 O 



(C27) 



nHlSfcTmH/vr) 1 ^ 2 ^ 



(C28) 



{Q)ti (E 1 / 2 )^, and (E 3 ^)^ were calculated for small graphite and silicate spheres for the 
interstellar radiation spectrum of Mezger, Mathis, & Panagia (1982) and Mathis, Mezger, 
& Panagia (1983). The results are presented in Table |] 
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Table 1: Idealized phases for interstellar matter. 



phase 


DC 


MC 


CNM 


WNM 


WIM 


%( cm~ 3 ) 


10 4 


300. 


30 


0.4 


0.1 


T(K) 


10. 


20. 


100. 


6000. 


8000. 


X 


io- 4 


0.01 


1. 


1. 


1. 


x n = n(H+)/n H 








0.0012 


0.1 


0.99 


x M = n(M + )/n H 


10~ 6 


0.0001 


0.0003 


0.0003 


0.001 


y = 2n(H 2 )/n H 


0.999 


0.99 


0. 


0. 


0. 



Table 2: Electric Dipole Moments for Selected Bonds' 



bond 


li 


aliphatic C-H 


0.3 


aliphatic C=0 


2.4 


aromatic C=0 


2.65 


aromatic C-OH 


1.6 


aliphatic C-OH 


1.7 


aromatic C-CH 3 


0.37 


aromatic C-C=CH 


0.7 


aliphatic C-C=CH 


0.9 


aromatic C-CHO 


2.96 


aliphatic C-CHO 


2.49 



a Dean (1992) 
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Table 3: Electric Dipole Moments for Selected Molecules 



molecule 


N 




P 


CsH 8 1 pentyne (trans) a 


13 


0.84 


0.23 


C7H5N benzonitrile a 


13 


4.18 


1.16 


HCnN cyanopolyyne fe 


13 


5.47 


1.52 


C 5 H 10 O cyclopentanoP 


16 


1.72 


0.43 


CgHio ethylbenzene a 


18 


0.59 


0.14 


CsHio orthoxylene a 


18 


0.62 


0.15 


Ci 2 H 10 acenaphthene a 


22 


0.85 


0.18 



a CRC Handbook (1997) 
b Bell et al. (1997) 



Table 4: Dust Grain Properties for Interstellar Starlight 



quantity 


graphite 


silicate 


cn-SiC 


(Q)r/a-r 


3.12 x 10~ 3 


5.79 x 10~ 4 


9.06 x 10~ 4 


(<5)W°-7 


8.92 x 10" 3 


3.20 x 10~ 3 


8.45 x 10~ 3 




2.56 eV 


4.63 eV 


9.41 eV 




2.78 eV 


5.28 eV 


9.68 eV 




3.03 eV 


5.86 eV 


9.87 eV 


g(100/im)/a_ 7 


1.93 x 10~ 3 


1.44 x 10~ 5 


9.03 x 10~ 7 


T d (ISRF) 


19.7 K 


17.4 K 


32.5 K 



Assuming interstellar starlight (Mathis, Mezger, & Panagia 1982; Mezger, Mathis & Panagia 



1983) with T 
See eq.( |C20HC23| ) 



1.525 xlO 10 cnrV 



UlSRF = 8.626 x 10" 13 ergs cm" 3 , and (E) T = 1.058 eV . 
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Fig. 1. — Charge distribution functions for grains of radii a = 5,7, 10, 15, 20 A for "Molecular 
Cloud" (MC), "Cold Neutral Medium" (CNM), "Warm Neutral Medium" (WNM), and 
"Warm Ionized Medium" (WIM) conditions. 
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Fig. 2. — Mode Z m , centroid (Z), rms charge (Z 2 ) 1 / 2 , and standard deviation o z for grains 
of radius a for CNM and WNM conditions. 
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Fig. 3. — Characteristic time scale tz [see eq. @] for changes in the grain charge Ze. Also 
shown is the characteristic rotational damping time tj [see eq. fl59|) 1 for a grain with charge 
Z m e. It is apparent that the approximation r z tj is excellent for all except the smallest 
(a < 4 A) grains. 



-44- 



10 s 



10 



0.1 



0.01 



10" 



i i i i i r 



i i r 



MC 




J I I I I I L 



io- 



io- 



a (cm) 



Fig. 4. — Dimensionless rotational drag functions for neutral collisions [F n , eq. (|TD|)], ion 
collisions [Fi, eq.(p0|)1, plasma drag [F p , eq.(p5|)1, and infrared emission [Fir, eq.( ^C| - |33| )1, for 
"Molecular Cloud" conditions. Plasma drag is more important than other gas collisional 
processes. Also shown is the ratio r H /r erf showing the relative importance of electric dipole 
damping. For molecular cloud conditions, electric dipole damping is not important. 
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Fig. 5.— Same as Fig. |4] but for "Cold Neutral Medium" conditions. Electric dipole 
damping dominates for a^4x 10~ 8 cm. Plasma drag F p is the dominant drag process for 
4 x 10~ 8 & a £ 1.5 x 10~ 7 cm. Infrared emission damping is dominant for a £ 1.5 x 10 -8 cm. 
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Fig. 6. — Same as Fig. § but for "Warm Neutral Medium" conditions. Electric dipole 
damping dominates for a ^ 1.5 x 10~ 8 cm, and IR damping dominates for larger grains. 
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Fig. 7. — Same as Fig. but for "Warm Ionized Medium" conditions. 
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Fig. 8. — Dimensionless rotational excitation functions for neutral collisions [G n , eq.(|38"D1, 
ion collisions [Gi, eq.(pP|)-(^TD], and plasma drag [G p , eq.([£J)] for CNM conditions. Ion 
collisions and plasma drag dominate, despite the low assumed ionization fraction of 10 -4 . 
Also shown is the "superthermal" excitation term G s for grains with recombination efficiency 
7 = 0.1 and N r given by eq.|54]. Because of the small assumed H fraction (1%) H2 formation 
torques are not important in molecular clouds. 
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Fig. 9. — Same as Fig. [8] but for CNM conditions. The contribution of infrared emission 
[(Sir, eq.(|44l)-(f46|)1 is also shown. Gm is the contribution if H 2 formation occurs randomly, 
rather than at preferred sites. G pe is the contribution from photoelectric emission (eq. |47D . 
Ion collisions and plasma drag dominate for a £ 10~ 7 cm, and systematic torques from H 2 
formation dominate for larger grains for the assumed H 2 formation efficiency 7 = 0.1. 
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Fig. 10. — Same as Fig. I and | but for WNM condit 10ns. 
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Fig. 11. — Same as Fig. | and | but for WNM condit ions. 
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Fig. 12. — Effective rotation rate uj ra( i = (u; 4 ) 1 / 4 as a function of grain radius a, for various 
environmental conditions, and with and without H2 formation. Also shown are the thermal 
rotation rates [eq. fljjD] at T = 20K and 8000K. The number of atoms N in a grain is 
indicated at the top of the figure. Also shown is the rotation rate for which the tensile stress 
equals 10 9 dyncm -2 . 
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Fig. 13. — Mean square 5u due to impacting ions relative to mean square u as a function of 
grain radius a, under CNM and WNM conditions. It is seen that the rotational excitation 
process may be treated as continous for a £ 7 A, but for smaller grains a single ion impact 
can result in an angular momentum much larger than the time average. For these small 
grains the angular velocity distribution will depart significantly from a Maxwellian. 
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Fig. 14. — Adopted sized distribution for small grains in diffuse regions (CNM,WNM,WIM), 
with 5% of the total carbon abundance in a log-normal component (see eq. |63|). In molecular 
gas, we assume the number of a ^ 3 x 1CT 7 cm grains to be reduced by a factor 5 (see text), 
so that the log-normal component contains 1% of the carbon. 
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Fig. 15. — Emissivity per H due to ultrasmall spinning dust grains under CNM and WNM 
conditions, for the size distribution of Fig. |14[ Also shown are observed emissivities from 
COBE D MR (open diamonds; Kogut et al. 1996); Saskatoon (filled diamonds; de Oliveira- 
Costa et al. 1997); and OVRO (filled squares; Leitch et al. 1997). 
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Fig. 16. — Same as Fig. [15] but for intrinsic electric dipole moments /ij reduced by a factor 
of two (/3b = 0.2debye). 
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Fig. 17. — Emissivity of CNM grains for 7 = and 7 = 0.1, where 7 is the fraction of 
incident H recombining to form H 2 . The H 2 formation torques have a minimal effect on the 
emissivity. 
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Fig. 18. — The central surface brightness of dust emission from a model cloud resembling 
L1755 (see text). Dashed lines show antenna temperatures of 100/xK and O.lmK . 
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Fig. 19. — The ratios a x /a and £ as a function of a. Here a is the radius of an equal- volume 
sphere, a x is the effective size from eq.(|3], and a is the factor by which the moment of inertia 
exceeds that of an equal- volume sphere. Tickmarks along the top of the figure show the 
number of atoms per grain. 



